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ABSTRACT : This report describes a method of using a digital computer 

root-locus program as a means for generating design charts which may 
be used to synthesize distributed- lumped-active networks. An appli- 
cation of the program is made to a specific network configuration in 
which the positions of the dominant complex conjugate poles and zeros 
may be chosen by the user. Design charts for this network are included. 
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AUTOMATED GENERATION OF DISTRIBUTED-LUMPED- 
ACTIVE NETWORK DESIGN CHARTS BY DIGITAL 
COMPUTER ROOT-LOCUS TECHNIQUES 

I. Introduction 

This is one of a series of reports describing the use of digital 
computational techniques in the analysis and synthesis of DLA (Distri- 
buted-Lumped-Active) networks. This class of networks consists of 
three distinct types of elements, namely, distributed elements (modeled 
by partial differential equations), lumped elements (modeled by alge- 
braic relations and ordinary differential equations), and active 
elements (modeled by algebraic relations). Such a characterization is 
applicable to a broad class of circuits, especially including those 
usually referred to as linear integrated circuits, since the fabrica- 
tion techniques for such circuits readily produce elements which may 
be modeled as distributed and active, as well as ones which may be con- 
sidered as lumped. The network functions which describe such circuits 
involve hyperbolic irrational functions of the complex frequency vari- 
able. For example, for the simple uniform distributed RC element shown 
in Fig. 1, the open-circuit voltage transfer function is : 



^1 cosh Jr C p 
o o K 



where p is the complex frequency variable and and are respectively 
the total resistance and the total capacitance of the element. The 
transcendental nature of the network functions of the type shown in (1) 




encountered in networks containing distributed elements produces many 
problems when the methods normally used in the analysis and synthesis 
of purely lumped (or lumped and active) networks are applied to the 
DLA network. The analysis problems may, in general, be overcome by 
using digital computational techniques to model the distributed ele- 
ments and also to solve the equations of the resulting network. A 
program for implementing such an analysis operation has been described 
in a previous report.^ - With respect to synthesis, one of the most 
fruitful approaches which has yet been proposed is the development of 
design charts which give the values for the network elements which will 
realize specific equivalent dominant poles and zeros. The points 
needed to construct these design charts may be found by matching the 
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characteristics of the networks at selected points along the jw axis, 

or by maximizing the magnitude of a network function at a given value 
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of complex frequency using optimization strategies. In this report a 
novel technique is developed for the direct production of such design 
charts. The method does not require the use of matching techniques or 
of optimization strategies. It uses a digital computational method for 
generating a root locus. When a set of such root loci are drawn on a 
single graph, the result is a complete design chart. The technique is 



quite general and it may be applied to a broad class of DLA networks. 

The report consists of two major sections. In sections II and III a 
discussion is given of the digital computer root-locus program devel- 
oped for this application. In section IV, the program is used to 
develop design charts for some specific DLA network configurations. 

II. Some General Aspects of a Root-Locus Program 

Several approaches to the use of digital computation techniques 
to plot root loci have been reported in the literature. Two of the 
more representative ones are described in the articles listed in the 
references 4,5. The root-locus program called ROOTLOC, which will be 
described in this and the following section of this report, has several 
features in which it differs from some or all of the root-loci deter- 
mining programs previously reported. Its major characteristics are: 

(1) A relatively simple algorithm is used to generate each of the root 
loci. Thus, it has been possible to keep the over-all size of the pro- 
gram quite small and, as a result, it is a relatively efficient program 
to compile and to execute. (2) A feature has been included which pro- 
vides for automatic internal scaling so that all of the points which 
determine the root locus are very nearly equispaced. (3) The function 
whose locus is to be plotted is described by a separate subprogram 
which is provided by the user. Thus it is possible to plot loci for 
transcendental functions as well as for the more usual polynomial functions. 

The digital computer program ROOTLOC is designed to construct a 
root locus plot of the function F(p,R) where p is the complex frequency 
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variable and R is the parameter which is being varied. Thus, for a 
given sequence of values of R, it is desired to find the corresponding 
values of p which satisfy the equation 

F(p,R) * 0. (2) 

The program requires that a user supplied subroutine called EQN be used 
to determine the value of the function of F(p,R) and the derivative 
function F’(p,R) where 


r (p.w - «) 

The basic algorithm used to follow the root locus may be expressed by 
assuming that some value of p, namely p^, is a point in the vicinity 
of a zero of F(p,R). If we now let be some new point in the vicin- 
ity of then the relation between an< * F(p^,R) may ke 

expressed by the relation 

F(p i+1’ R) = F( V R) + <P i+1 - P ± > ^'(Pi.R) + ••• (4) 

If the new point p^ + ^ ^ as ^ een chosen so that it is actually a zero of 
F(p,R), then the left member of the above equation is zero and, ignor- 
ing terms of higher than the first order in the right member, we may write 

F(p 1 ,R) 

p i+l " P i " F 7 (p i , R) 


( 5 ) 
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The relation given above may be applied iteratively as an algorithm to 
locate a zero of F(p,R) starting from some point which is close to such 
a zero. This equation is, of course, the well-known Newton-Raphson 
relationship. A discussion of its convergence properties may be found 
in any of the standard texts of numerical analysis. The problem 
solved by this algorithm is, of course, exactly the problem which is 
encountered when a point on a root locus for a given value of R is 
known and it is desired to find another nearby point, also on the root 
locus, representing, a slightly changed value of R. Thus the relation 
of (5) is directly applicable to the root-locus problem being considered 
here. 

The general operation of the ROOTLOC program is readily described. 
The user supplies the program with a known starting point, i.e., a 
value of p which is a root of F(p,R) • 0. The program then makes a 
change in the parameter R and iteratively applies (5) until convergence 
is found, i.e., until a new value of p is found which satisfies (2). 

This new value of p is thus a point on the root locus. Next a check is 
made to insure that the point which has been found is a convenient dis- 
tance (for plotting) away from the preceding point. If it is not, the 
program automatically makes an appropriate change in the value of R and 
the Newton-Raphson algorithm is again applied. This procedure ensures 
that the points defining the root locus will be equispaced, even though 
the sensitivities of the roots with respect to R may vary widely at 
different parts of the root locus. The process is continued until some 
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terminal value of R, specified by the user, is reached. As many sets 
of starting points for the root loci, and as many ranges of the para- 
meter R may be used as are desired. Thus, different branches of the 
overall root-locus plot are easily found by repeated application of the 
basic logic described above. After all the points determining the root 
loci! have been located, they are plotted using an internal plotting 
capability which is included as part of the program. This plotting 
subroutine assigns alphabetical symbols to the various points as they 
are computed. These alphabetical symbols are then reproduced directly 
on the plot so that the value of the parameter R corresponding with any 
point on the root locus may be readily determined from an accompanying 
tabulation of data. Multiple use of each alphabetical character is 
made if the number of points which is plotted exceeds 25 (the character 
I is not used). Thus, if a hundred points are to be plotted, the first 
4 will be plotted using the character A, the next 4 the character B, 
etc. Since the relations for F (p , R) are supplied by the user, the pro- 
gram may be used to plot the loci of transcendental functions of the 
kind encountered in DLA networks as well as the more conventional poly- 
nomial functions. A simplified flow chart of the sequence of operations 
used by the ROOTLOC program is given in Fig. 2. Some more detailed 
information on the operation of the ROOTLOC program is given in the 
next section. 

III. Operation of the Dig! tal Computer Projtram ROOTLOC 

The digital computer root-locus plotting program ROOTLOC is 
designed to implement the general computational process described in 
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the preceding section. The program requires the use of two subroutines. 
The first of these is a plotting subroutine entitled XZPLT. This is 
supplied as part of the root- locus program package. The second subrou- 
tine which is called EQN, is provided by the user. The details of 
these two subroutines follow: 

XZPLT — Subroutine for making an x-y plot of a set of data points with 
"z" coordinate information specified through the use of various 
alphabetic symbols to indicate the values of the parameter R. The 
subroutine plots the abscissa vertically downward on the printed 
page. Thus, as large a total range of abscissa scale as is desired 
may be used. The scale used for the abscissa is 10 units per inch. 
Since there are six printed lines per inch on a standard digital 
computer printer, a range of 90 units will cover 55 printed lines. 
This corresponds approximately with a single printed page. The 
ordinate is printed across the page. Again, a scale of 10 units 
per inch is used. Since there are 10 characters per inch on the 
standard printer, a range of 100 units on the ordinate scale encom- 
passes 100 characters of printing and provides a 10 inch wide plot. 
A listing of the subroutine XZPLT is included with the listing of 
the main program given in Appendix A. A summary of its character- 
istics and a basic description of its operation can be found in 
the literature.^ 

EQN — A user supplied subroutine which computes the value of function 
F(p,R) and the value of the derivative function F'(p,R). The 
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argument listing of this subroutine is (R, P, F, FP), where R is 
the (real) value of the parameter being varied, P is the (complex) 
value of the complex frequency variable, F is the (complex) value 
of F(p,R) and FP is the (complex) value of F’(p,R)« Thus R and P 
are inputs, and F and FP are outputs. 

The main program for ROOTLOG has three major functions. These are: 
(1) To read the input data; (2) to compute a succession of points which 
define the root locus; and (3) to plot these points on a two-dimensional 
graph. The details concerning these three functions are covered in the 
remainder of this section. 

The first function of the digital computer program ROOTLOC is to 
read the pertinent data indicated by the program. There are three types 
of data. The first type is a title statement which identifies the pro- 
gram. This is punched on a single card in alphanumeric format and may 
be from 0-80 characters in length. This title information is only read 
once for any given submission of the program, and it is reproduced on 
all the program output. The second type of data which is read by the 
program is the values of the 10 internal parameters which are used in 
the program. These parameters are read only once. They are identified 
by their FORTRAN variable name. Their function is explained in the list 
which follows: 

FMAX - This variable is a convergence factor which is used as a 
criterion to decide when the Newton-Raphson algorithm has found a 
value of p sufficiently close to an actual zero of F(p,R), It is 
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the maximum value which the quantity F(p,R) is allowed to have for 
any given set of values of p and R. The program will continue to 
apply the basic relation given in (5) until the absolute value of 
F(p,R) is less than this value or until the maximum permitted num- 
ber of iterations (see ITMAX below) is exceeded. A value for this 
parameter which has been found suitable for many applications of 
the program is .001. 

SCALE - This parameter is a scaling factor which is applied to the 
locations of the data points as determined by the R00TL0C program 
before they are plotted. This factor must be specified by the 
user since the plotting subroutine XZPLT automatically assigns a 
range of 100 units for the "y" coordinate. For example, if it is 
desired to have the ordinate range start at the "x" axis, then the 
maximum ordinate value permitted is 100. The scaling factor must 
be chosen so that the data points will lie within the appropriate 
hundred-unit range which is chosen for the ordinate. Similarly, 
the scaling factor must be correlated to the "x" axis values which 
are to be plotted. The total range of the abscissa or "x" axis is 
completely variable* However, as pointed out above, a range of 
approximately 90 units requires one page of printed output. Thus, 
if a one-page plot is desired the scaling should be so chosen as 
to insure that the scaled data points fall within a range of 90 
units for the abscissa. 

ITMAX - This parameter specifies the maximum number of iterations 
which is permitted for the Newton-Raphson algorithm. If this 
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number of iterations is exceeded without convergence (as specified 
by FMAX) having been obtained, the program will stop computation 
of the root locus which it is currently following, and will go on 
to other root loci. If no other root loci are to be computed, it 
will proceed to the plotting phase of the program. A value of this 
parameter which has been found suitable for many of the applica- 
tions for which this program has been used is 10. 

IRMAX - This parameter limits the maximum number of iterations 
which are permitted in the algorithm which adjusts the change made 
in the parameter R so that the computed data points are spaced 
equally in the resulting plot. If the number of adjustments exceeds 
this value the program will stop its efforts to compute the root 
locus which it is currently following and will go on to other root 
loci, or, if no other root loci are to be computed, it will proceed 
to the plotting phase of the program. A value of this parameter 
which has been found suitable for many of the applications for 
which this program has been used is 10. 

NSX - This parameter specifies the maximum abscissa value to be 
used on the plot produced by the XZPLT subroutine. 

NSY - This parameter specifies the maximum ordinate value which is 
desired on the plot produced by the XZPLT subroutine. The minimum 
value which will be used on the plot is automatically set 100 units 
lower than the value specified for NSY. 
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NNP - This parameter specifies the range of abscissa values which 
it is desired to have on the plot. If this parameter is set to a 
value of 90 then the plot will cover approximately one page. 

LMX - This parameter specifies the maximum number of points that 
it is desired to plot. The maximum value for this parameter which 
is allowed by the dimensioning of the program is 200. 

LPLT - This parameter provides an option by means of which the 
plotting capability of the R00TL0C program may be suppressed. 
Normally it is set to 0 for the case in which a plot is desired. 
However, if the parameter is set to 1, the plot will be suppressed. 
LPRNT - This parameter is used to select either of two formats for 
the printing of data as it is computed. If the parameter is set 
to 0 the results of each step in the various iterative processes 
of the program will be printed. If the parameter is set to 1, a 
reduced output is used in which only the final (successful) results 
of the iterative processes are printed. 

The third type of data which is read by the R00TL0C program is the 
(complex) value of the starting point which is chosen for the root of 
the F(p,R), the initial and final values of the parameter R, and an 
initial value for the variable DR which is used in changing the para- 
meter R. As many sets of this third type of data may be used as is 
desired. Each such set will be used to construct one root locus, or 
one section of a root locus. This makes it possible to avoid computa- 
tions at branch points, where F'(p,R) is zero, and thus where the 
iterative relation given in (5) will not converge. 
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A detailed set of instructions for preparing input data for the 
digital computer program ROOTLOC is given in the Appendix, together 
with a flow chart and a listing of the program. 

IV. Application of the Program ROOTLOC to some DLA Networks. 

As a verification of the utility of the approach implemented by 
the ROOTLOC program described in the preceding sections of this report, 
the program has been applied to analyze the well known low-pass DLA net- 
work configuration shown in Fig. 3. The voltage transfer function for 
2 

this network is 


V I v ' 

1 1 - cosh© + — cosh© 

where 0 ■ yRC p, and where R and C are respectively the total resis- 

oo r o o v J 

tance and total capacitance of the uniform distributed RC line, K is 
the gain of the VCVS represented by the triangle shown in the figure, 
and p is the complex frequency variable. 



Fig. 3 
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If the basic fundamentals of root locus techniques are applied to 

the function given above (these techniques as conventionally defined 
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must be extended so that they apply to transcendental functions ) for 
an R q C o product of 10, the resulting (hand) computations produce the 
root- locus shown in Fig. 4 in which the movement of both the dominant 
and non-dominant poles of the network function are illustrated. The 
arrows shown in the figure indicate the direction of increasing K. 

In order to apply R00TL0C to develop a plot of the type shown in 
Fig. 4, a subroutine EQN and a set of data cards are required. List- 
ings of these are shown in Fig, 5 and Fig. 6 respectively. The output 
data from the R00TL0C program giving the scaled values of the points 
on the root locus, the values of gain and the associated alphabetic 
symbol used for plotting are shown in Fig. 7. The resulting plot is 
shown in Fig. 8. It is readily verified that this plot agrees with the 
one shown in Fig. 4. 

In the above example, we have verified the feasibility of applying 
the R00TL0C program to determine a root locus for a DLA network config- 
uration. Now let us see how this program may be applied as a synthesis 
tool to develop design charts for specific DLA network configurations. 

As an example of such an application, consider the network shown in 
Fig. 9. It is well known that this network will realize a voltage trans- 
fer function of the form 



2 2 

P + b-jP + b 2 
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Figure 5 Subroutine containing equations for the example problem 
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over certain ranges of the constants a and b 2> We will now show that 
this circuit can also be used to synthesize transfer functions of the 
form^ 


V 2 G(p^ + a^p + a 2 ) 

BS ' 1 1 1 

V 2 

V I p Z + b x p + b 2 


( 8 ) 


To see this let us first consider the voltage transfer function for the 
network. This is 


V, r c Q sinh© + R 

2 f o 

V 1 R - R cosh© + \ , . . °l R o s 0sinh0 

o o K { (r f + — — ) 

o 

(?) 

+ (R + 2C_r c R p) cosh© - 2C T r c R p} 
o I f o r I f o 


where 

0 - 7r C p 
v o o 


( 10 ) 
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The numerator of this voltage transfer function contains three varia- 
bles, R , C , and r_. If R and C are chosen such that a zero is 
’ o’ o’ f o o 

possible for a=0, w=l, then, as r^ is varied the zeros of the voltage 
transfer function can be located off the jw axis in the right or left 
half of the complex plane by varying the feedback resistance r^ as 
shown in Fig. 10. 

After the zeros of the voltage transfer function have been chosen 
the values of R q , C o » and r^ are fixed. Looking at the denominator of 
(9) it can be seen that this leaves two variables in each circuit that 
are still free to be chosen. These two variables are gain (K) and the 
value of the input capacitor (C^) . Thus, for fixed values of either of 
these variables we can generate root loci as a function of the other 
variable for the specified values of R q , C , and r^. Listings of the 
subroutine EQN for the two cases are shown in Figs. 11 and 12. Fig. 11 
covers the case where r^ is the parameter and Fig. 12 covers the case 
where K is the parameter. Some root locus plots for various cases have 
been superimposed in Figs. 13, 14, and 15. These figures provide a set 
of design charts which are applicable to the indicated regions of the 
complex frequency plane. 

The use of the capacitative load at the input to the VCVS as shown 
in Fig. 9 does not permit the poles of the resulting transfer function 
to be located in all regions of the complex frequency plane. Additional 
coverage is provided by replacing the shunt capacitor with an input 
resistor R^. The voltage transfer function for this network is 
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RETURN 
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Figure 11 Subroutine containing the equations of the characteristic equation 


of the open circuit transfer function- of the circuit shown in Figured 


when r^ is the varying parameter. 
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—Filiure 12 Subroutine containing the equations of the characteristic equation 


of the open circuit transfer function of tha circuit shown- in Figure, 9 

when the gain-(K) is -the varying parameter. 










Figure 15 Family of root loci for the characteristic equation of 
the open circuit voltage transfer function of the above SLA network 
as K and vary. 
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r.0 sinhG + R 
f o 


R - R coshG + ~ { (r. + R o 
o o K f — 


Vo" 


) 0 sinhG 


2R r- 2R r, 

+ (R + ) coshG - } 


R-r 




( 11 ) 


Note that this function has the same numerator as the one given in (9). 

Thus, if we choose values of R , C , and r, such that the zeros of the 

o o r 

voltage transfer function are specified, there will again be two unspec- 
ified variables which can be used to determine the pole locations. 

These variables are the gain K and the input resistance R^. Thus, for 
fixed values of either variable we can generate root loci which are 
functions of the other variable. Listings of a subroutine EQN which 
will implement the two cases are shown in Figs. 16 and 17. Some design 
charts which have been prepared by superimposing several root-locus 
configurations for this network are shown in Figs. 18 and 19. If we 
compare the various design charts which have been presented so far, we 
see that, in general, the shunt input capacitor network shown in Fig. 9 
is suitable for network functions in which the poles are closer to the 
origin than the zeros, while the shunt input resistor network is suita- 
ble for network functions in which the poles are farther from the origin 
than the zeros. 

In this section we have illustrated applications of the ROOTLOC 
program to develop design charts for some well known DLA network con- 
figurations in such a way as to extend the capabilities of these 
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Figure _16_ Subroutine containing^ the equations of the characteristic equation 

g£ the open circuit transfer ftoction of the circuit shown in Figure! 9 ; 

•when rf is the varying parameter. (Input C replaced by input R) 
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configurations from the original case in which the zeros of the voltage 
transfer function were permitted only on the ju) axis to the case where 
zeros are permitted anywhere in the complex frequency plane. An exam- 
ple of the synthesis of a special network function using these design 
charts is given in the following section. 

V. Example of Design Procedure 

As a general procedure for the use of the design charts developed 
in the preceding section, to synthesize a network function having a 
dominant pair of complex conjugate poles and zeros, the following steps 
may be followed: 

1. Find the zero on the normalized numerator root locus shown in 
Fig. 10 with the same damping ratio as the complex zeros of 
the filter being synthesized. This determines the value of r^ 
for the normalized filter and the frequency normalization fac- 
tor (a , where 

n 

real component of zero of filter 
w . being synthesized 
n real component of zero on normalized 
numerator root locus with same 
damping ratio 

2. Frequency normalize the desired pole locations for the network 
being synthesized using the factor 

3. Check the families of curves shown in Figs. 13, 14, 15, 18, 
and 19 at the desired normalized pole location to find the 
one with the desired r^ curve passing through this normalized 
pole location. This determines K and Rj. (or C^.) for the 
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frequency-normalized filter. If the desired r^ curve does not 
pass through the desired pole location on any of these families 
of curves, then interpolation between families of curves will 
be necessary. 

4. The values of the components for the filter being synthesized 
are found from the values of the components for the frequency 
normalized filter by dividing the values of the frequency 
dependent components by 

As an example of the application of this procedure, consider the 
voltage transfer function 

V 2 (p) (p + ,3 + U.5)(p + .3 - il.5) p 2 + ,6p + 2.34 

V p> ‘ (p + 1.6 + J1.6) (p + 1.6 - 11.6) - ^2 + 2 + 512 

Applying the steps defined above to this network function we obtain 
the following results from each of the steps. 

1. The feedback resistor r^ ■ 1,9. The frequency normalization 

factor to = 2. 
n 

2. The normalized roots of the denominator are p = -0.8+j0.8. 

3. The family of curves for which the point (-0.8 + j0.8) is on 
the r^ = 1.9 curve is Rj. =■ 10, and the gain of the VGVS is 1.3. 

4. The DLA network for the desired filter is shown in Fig. 20. 

It should be noted that, in the example given above, an exact solution 
was found on one of the design charts developed in the preceding section. 



33 



This, of course, will not always be true. For a more general synthesis 
capability it would be desirable to generate additional families of 
curves. However, the curves given in this report can be used as the 
basis for interpolation to approximate the values of the network ele- 
ments which will produce a broad range of the most commonly encountered 
network functions. 
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APPENDIX 


Instructions for Preparing Input Data for ROOTLOC Program 


Card 


No. 

Variable 

Columns 

Format 

Purpose 

i 

LTR 

1-80 

8A10 

Title card for problem 

2 

FMAX 

1-10 

E10.0 

Convergence factor for 
Newton-Raphson algorithm 


SCALE 

11-20 

E10.0 

Scaling factor for locations 
of roots of function 


ITMAX 

21-23 

13 

Maximum number of iterations 
permitted for Newton-Raphson 
algorithm 


IRMAX 

24-26 

13 

Maximum number of iterations 
permitted for algorithm 
which adjusts change in 
parameter R 


NSX 

27-29 

13 

Maximum abscissa value 
desired on plot 


NSY 

30-32 

13 

Maximum ordinate value 
desired on plot (minimum 
value is 100 units lower) 


NNP 

33-35 

13 

Range of abscissa values 
desired on plot (90 covers 
one page) 


LMX 

36-38 

13 

Maximum number of points 
that it is desired, to plot 
(if not specified, set to 
200 by program) 


LPLT 

39-41 

. 13 

Indicator for plotting: 0 
if plot desired; 1 if no 
plot desired 
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lprnt 

42-44 

13 

Indicator for printing data: 
0 if data is desired for 
each step of process; 1 if 
only successful steps are to 
be printed 

3 

P (Real) 

1-10 

ElO.O 

Starting value of 


P(Imag) 

11-20 

E10.0 

root location 


RA 

21-30 

ElO.O 

Initial value of parameter 
R which is to be varied 


RB 

31-40 

ElO.O 

Final value of parameter 
R which is to be varied 


DR 

41-50 

ElO.O 

Initial change to be tried 
in changing the parameter R 


fmax 

51-60 

ElO.O 

New value for convergence 
factor (if no new value is 
specified here, the original 
value read on card no. 2 is 
retained) 

4 

Blank card i 

signifying 

end of data 

Rotes : 

1. As 

many cards of 

type No. I 

3 may b e 

used as is desired before the 


blank card. The program may thus be used to plot loci starting from 
different root locations and/or for different ranges of the parameter 
R. If the number of points computed exceeds LMAX (maximum value 200), 
the program will ignore any remaining cards and proceed with its plot- 
ting cycle. 

2. If it is desired to have the user-supplied subroutine EQN read 
input data, the subroutine must be arranged so that it does this only 
on its first call. The data cards should be inserted following the 
first Card No. 3 which is used in the data card deck. 



Read FMAX, 


Start 

i 

Read and print LTR 

4 / 

Set L » 1 

A 

SCALE, ITMAX, IRMAX, LMX, 

4 ' 


LPLT, LPRNT 


Is LMAX 0? 


No 



Yes 


> LMAX = 200 



Print RA, RB, FMAX 


i 

Set R m RA 

> i 

V 



Call EQN to find F and FP 
Set FM = |F[ and FPM - | FP | 


X(L) * Re (P) x SCALE 
Y(L) = IK (P) x SCALE 
Z(L) = R 



Flow Chart for Root Locus Program 



Is LPRNT = 0? 



Print L, IT, IR,P,F,FM,FP,FPM, R,Dr,DZ 

-> l 

Is FM > FMAX? 


Is INR > 0? 


IT - 0 

i 

Is R = RA?- 


Is R * RB? 


Is DZ > 6? 


Is DZ > 4? 


DZ=l-f-- Is DZ < .5? 

i Yes | 


Is DZ < 2? 


IR - 0 

4 ' 

Is LPRNT m Oi - 


ls FPM < .001? > Multiply 

FP by 

No 100 


P * P - F/FP 

4 

INR = 0 

4 

IT - IT + 1 

4 

Is IT ^ ITMAX 


Print L, IT,IR,P,F,FM,FP,FPM, R,DR,DZ 

© , ttr 

v — ^ L «■ L + 1 


Flow Chart for Root Locus Program (page 2) 




Is LPLT - I? 


Call XZPLT 


DR - change of variable 

DZ - distance from preceding point 

F - value of function 

FMAX - convergence factor for Newton -Raphson algorithm 
FP - value of derivative of function 

INR - indicator to insure use of Newton -Raphson after each gain change 

IR - counter for gain adjustments 

IRMAX - maximum number of gain adjustment cycles 

IT - counter for Newton -Raphson cycles 

UMAX - maximum number of Newton -Raphson cycles 

L - counter for point being computed 

LMX - maximum number of points to be computed 

LPLT - plotting indicator (1 if no plot desired) 

LPRNT - reduced printing indicator (1 if reduced printing desired) 
LTR(I) - title 
P - position of root 
R - variable 

RA - initial value of variable 
RB - final value of variable 
SCALE - scale for plotting 


< RB?— 

Yes 

\ 

Is 

No 

( 

» RB 

1 

— — 

R 





500003 
>00003 
0 0 3 
| vu 0 1 } 

|ooon 

00001 7 
•000017 
Jqoopq 
Iooo2i 

OOOgBl 

boooS3 
1)0 0053 
000071 


000071 

000075 


000075_ 
000076 
b 0 0 0 7 6 
boo i oo 
000115 
000116 
pooi?o„ 
000120 
000121 
p^l33~ 

000133" 
1)00135 
boo 140 
000142 
000144 

booiso 

000152" 
000154_ 
000155 
000157 
000165' 
|0 0 0 1 7 0 
1000222" 
0 00222_ 
|0 0 0226 
1000230 
*0 0 02 3 0 
000232 
|o 0 0 2 3 4 
1000240 
0 00 244 
1000247 
1000252' 
000253 
"1 254* 

|o u 0 3 0 6 


PROGRAM MAIN ( INPUT ,001 PUl ) 

MAIN PR0OHAM f-oR HOU I LOCHS PLOf 
COriPLF X P , F , F P 

DIMENSION LTR(o),X(2u0)»Y(200) »Z (200) 

READ 1 05* UH 

1 05 FORMAT (HA 10) 40 

PRINT 110, LIP ' ; \\ 

llO FORMAT (1H1,8A10) 

UZ=:0. 

L s 1 

RF.AO 1 15»KMAX,bCAlt» UMAX* IRMAX»NSX,NSY»NNP»LMX»LRLT.LPRNT 
IF (LMX.LU.O) L M X = 2 0 0 
115 FORMAT <2E10.0*HI3) 

PRINT 1?0»F'MAX»1TMAX, I RMAX . LMX , 5C ALE 
1 20 FORMAT ( 1. HO » 22HM { N VALUt OF F UNC UON* ♦ E9 . 2 , 5X ♦ 

VibHMAX I iERAl 1UNS«*I3*5X* 16HMAX AU JUSTMtN I S* » I 3 
1 * 5X » 1 ) HMAX POINTS** I3*5X,6HSCALt = »E8.1//) 

PRINT 125 

125 FORMAT ( 1H0*2X»1HL» 1X*?HIT» 1X»2H1Ri9X*4HR00T* 15X»8HFUNCTI0N*8X* 

1 5 H ( M A G ) |9X,10HUER1VATIVE*8X,5H(MA6) * 6X * 1HR »9X 1 2H0R *8X , 2HUZ/ ) 

__130_ 1« = 0 

I N R s 0 
I TsO 

READ 135»P*RA*RB«0R*FZ 
D Z e 0 * 

"“"IF (FZ.NE.O.) FMAX*FZ 

135 FORMAT (BEiQ.O) _ . 

IF (OR.EQ.O.) <jO TO 180 
PRINT 1 40 *RA » Rb * FMAX 

14 OF ORMAT ( 1H0 *22flRANGE OF VARIABLE FROM , El 0 . 2 » 2X , 
_12HTO*E10.2*SX*22HMIN VALUE OF FUNCTION* *£9.2/) 

' R*R A 

142 C ALL EON (R«P*F»FP) 

FMsCABS ( F ) 

FPM*CAR$(FP) 

X(L)*REAL(P»*SCALE 

Y (L) sAIMAG(P) ASCALE 

Z(L)*R 

IF (R.EG.HA) GO TO )44 

IF (R.EQ.Rti) GU TO 144 

OZsSQRT ( (X (L)-X(L-l) > **2* < Y ( L > - Y ( L- 1 > > **2 ) 

IF (LPRN r .GT « 0) GO TO ) 46 
_1 44 PRINT 145*L* 1T»IR.P*F.FM,FP,FPM*R»0R»DZ 

~ '145 'FORMAT UX*3I3*2<E10.?*2H*J*E^#2) » 2EI 0 . 2 * 2H* J.E9, 2 , 4E 1 0 .2 ) 

146 IF (FM,GT,FMAX? GO TO 170 __ 

IF ( I MR , G T.O) GO TO 170 
I TsO 

IF (R.FQ.RA) GO TO 150 
IF TR.EO.Rb) GO TO 150 
IF" (OZ.G1 .6) l)Z*5.0 
IF (DZ.GT.4.) gO TO 160 

if (dz.lT . 0.5) oz«i.o - 

IF (DZ.LT. 2.) GO TO 160 
150 19*0 " 

IF (LPRNT.EQ.O) GO TO 152 ’ 

PRINT 145* L» IT »1R.P.F»FM»FP*FPM»R*DR»DZ 
152 L*L*l ‘ , 'CM 





. 

00031 U 


IF ( L * (IT . LMX ) 00 TU 1H0 

.. J 

0 0 0 J 1 3 


IF (W.Ftl.Hb) 00 TO 130 


00031b 


KS«*l)W 

j 

00031b 


1NMS1 

J 

000317 

1 5 1 

IF (OP,GF.O.) 00 10 Ibb 


003? 1 


IF (R.GT.Rh) GO TO 142 

"1 

0 0 0 3 ? 4 


H = R H 

1 

000324 


00 10 142 


0 0 03?b 

lbb 

IF (H.LT.RH) 00 TO 14? 

ii 

000330 


R=PR 

j 

0 0 u 3 3 0 


00 TO 142 


1 000331 

160 

R=R-OR 


000333 


i)RsJ,06*UH/OZ 

1 

00033b 


H*H+OR 


00033b 


l MRS 1 


000337 


iRsiWM 

— 

[ 000341" 


IF ( I P , L L • IRMA A ) GO TO 151 


| 000343 


GO TO 1 30 


00034* 

1 7 0 

IF' (FPM.LT.O.QUl ) FP*FPMlOO f 

71 

0003bb 


P=P~F/FP 

j 

000367 


1 T s 1 T ♦ 1 


000370 


IMP so 


0 00371 


IF ( IT.Ltt ITMAA) GO TO 142 


000373 


GO TO 130 


000374 

180 

L=L-1 


000376 


IF (LPLT »EQ. 1 ) STOP 


000401 


call XZPLT <L*A.Y*NSX» NSYiNNP 

a j 

000410 


STOP 


000412 


END 



) 0 0 0 1 0 
) 0 0 0 1 0 


h ouul o 
booo 1 4 
000u?0 

poo 0?0 

f)00u?4 
000024 
O0O02B 
f)0002b 
0 0 0 0 3 3 


m :• 

' 

■' 

&&S ■ 


000204 


bUbHoin JNE X/i’LT (ND.X. f * NisX * N5 Y * NNP » Z ) 

u 1 MENS I OH X ( 2 0 U ) » Y ( 2 0 0 ) « L d 1} * L I N't ( 10 1 ) * Z ( 2 0 0 > ♦ JL < 25 ) » M ( 2 0 0 ) 
uATA ( JL ( 1 ) » 1*1 *25) / IHA » lnw» 1 HC» 1 HU « lHE* IMF t 1HG* lhH» 1 HJ* lhK i 
liHL » 1 HM ♦ l H Si » 1HU* 1HP » 1 HQ * 1 HH * 1 Hb * i H f » IMUf 1 HV * 1 HW » 1 HA » 1 H Y ♦ 1 HZ / « 

2 J M / 1 H - / * J P / 1 H ♦ / ♦ J 1 / 1 H % / » Jtf L A NK/ 1 H ; / » Jl/l H 1 / 

N O A - ( MD ♦ 2 4 ) / 2 5 
2M PRINT ?2<1 

22 0 f-ORMAT ( IH1 , 10A» 19HDATA FOR X-Y-Z PLOT) 

ph (N r ?.?.b llllillllllliillillllllllllil^ 

225 F ORMAT ( IH0*15X« IHX« 19X « 1 HY , 19X ♦ 1 HZ » IbX » 6HSYMHOL/ ) 

L Za 1 
J=1 

227 UO 24 0 1 = 1 * ND A 

PRINT 235 1 J » X ( J ) ,Y < J) *7 ( J) » JL <LZ> 


0 0 o U 7 1 

235 

FORMAT ( 15 * 3 L 20 . 

000071 


m ( J ) = JL (LZ> 

0 0 0 0 7b 


jsj+ 1 

000076 


IF" ( J-MO) 240 * 240 * 

000104 

240 

CONTINUE 

000107 

245 

LZsLZD 

0001 n 


GO TO 227 

000111 

“250 

NOMsND -1 

000113 


UO 200 I» 1 *NQM 

000114 


I A* I ♦ 1 

000116 


UO 200 J*I AtND 

000120 


IF(X(I)-XU)> 200 * 

Jo 00 124 

215 

TEMPcX(I) 

000126 


A ( I ) ®X ( J ) 

1000132 


A ( J ) sTEMP 

|roi 33 


TEMP=Y ( I ) 

v. j\ 35 


Y(I)aY(J) 

,000140 


Y ( J) =TEMP 

S 000 14 1 


(DEMPsM ( 1 ) 

000 143 


M ( I ) sM ( J ) 

000146 


M ( J) sKTEMP 

000161 


Z < I ) =2 ( J ) 

000154 


2 ( J ) -T EMP 

000156 

200 

continue 

000163 


NPs (NNP/ 10) *6 

000 16 b 


' XNPSNP 

000170 


aNS 3 (NSX/.l 0 ) #6 

-000174 


YNS»NSY 

000176 


00 101 1 = 1*11 

J 00(>177 


Ld)al 0 *l- 110 >NSY 

000203 

101 

CONTINUE 


200*215 


PRINT 1 02 »L 


000212 

102 

FORMAT (lHl,2Xfll 

000212 


N = 0 

000213 


K=1 

000214 

105 

NOsO 

000215 


UO 110 1=1*10 

000222 


nOsNQ* 1 

000224 


LINE (NQ)SJP 

000226 


UO 110 J- 1 *9 

000230 


NQaNQ+1 

7232 

110 

LINE (NO) aJM 

0u0240 


LlNE(101)aJP 




llliillil! 



|000261 


IF ( N ) U 9 ♦ 11 5 , U 9 

000243 

l 1 s 

IvXrA (K ) 0.0-XNS + a.NP* ,499999 

000251 


IF ( N a ) 1)6*119*117 

000253 

1 1ft 

I'jAsQ 

000254 


30 T(.) 119 

Ocbb 

117 

IF < NX-NP ) 119*119,118 

j v 00260 

1 1 8 

NXsNR 

*000262 

119 

J, F (NX-N) 125*120*125 

000264 

1 2f) 

ijYsY (K) ♦ 101 ,499999-YNS 

10002/0 


IF (NY-1) 122,126*123 

j 0002 /3 

122 

L I ME. ( 1 ) =J2 

0002/5 


00 TO 12/ 

, 0002/5 

123 

IF (NY-101) 126*126,124 

1 00030U 

U4 

LINE ( 101 ) =J2 

*000302 


00 TO 12/ 

000302 

12ft 

LINE<NY)sM<K) 

[000306 

127 

iA*K * 1 

J 000310 


IF (K-NO) 115*115*128 

000312 

125 

IF (N) 128*135*128 

j 000313 

12R 

IF (M/ft- (N-l ) /ft) 130*130*135 

j .0 0 0 3 2 2 

* 1 3 0 

PRINT 131* L I Nt 

000330 

.131 

FORMAT (5X* 10 1 A 1 ) 

, 000330 


*60 TO 140 

I 000334 

135 

l\Nc (N»10) /6* NSX-NNP 

1 000342 


PRINT 132* NN*LINE 

000351 

13? 

FORMAT (1X.I4»101A1) 

000351 

140 

IF (N-NP) 145*160*160 

000357 

145 

N = N ♦ 1 

000361 


IF (N/6- (N-l) /6) 148*148*105 

| 000367 

148 

00 150 1=1*101 

, "40371 ' ' 


L I NE ( I ) a JBLANK 

.0373 

15 0 

CONTINUE 

, 000375 


DO 152 1 = 1*101*10 

i 000376 


L I NE ( I ) a J I 

J 000400 

15? 

continue 

000402 


GO TO 119 

[000402* " 

160 

RETURN 

J 000403 


END 



